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SUMMARY _ . - 


A number of calculations of bending- tor slon wing 
flutter are made at two Mach, numbers, M ~ 0 (incom- 
pressible case) and M = 0 . 7 , and the results are com- 
pared. The air forces employed for the case of M- = 0.7 
are based on Frazer’s recalculation of Possio's results, 
which are derived on the assumption of small disturbances 
to the main flow. For ordinary wings of normal density 
and of low bending frequency in comparison with torsion 
frequency, the compressibility correction to the flutter 
speed appears^to be of the order of a few percent; whereas 
the correction to the flutter speed for high-density wing 
sections, such as propeller sections, and to the wing- 
divergence speed In general may be based on a rule using 

the (l - M 2 ) 1 ^' factor and, for M = 0-7, represents 
a. decrease of the order of 17 percent. 


INTRODUCTION 


The question of the Influence of the compressible 
properties of a gas on wing flutter is,, of course, 
directly tied to the primary problem of determining the 
air forces and moments on oscillating airfoils moving at 
high forward speeds. 'This problem has been attacked by 
Possio (reference 1), along lines indicated by Prandtl, by 
a procedure utilizing the pressure or acceleration poten- 
tial and 'the method of linearization of the., equation sat- 
isfied by the acceleration potential for small disturb- 
ances 'to the main flow. A review and summary 'of Foasio’s 
work with certain simplifications have been given by 
Frazer (reference 2 ) . Frazer and Skan (reference 5) 
listed improved numerical tables of Possio’s results and 
made some numerical applications to the flutter problem. 
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It appeared worth while to perform additional cal- ’ 
culations along similar lines utilizing the notations and 
parameters more familiar in this country. The present 
report has for its limited objective the reporting of the 
results of a number of pertinent calculations on bending- 
torsion flutter for* a stream Mach number M = 0.7 and 
the comparison of - these results with those given in refer- 
ence I). for the incompressible case based on the Theodorsen 
theory (reference 5 )* 

The numerical accuracy of the results of— Poaslo's 
•theory and method deteriorates as M approaches unity 
and a3 the frequency increases. It -has been estimated 
that the theory is not safely-aoplicable much beyond M = Q-7 
nor at M = O.f for values of the reduced frequency k 
much beyond 1. Thus, the transonic or supersonic ranges 
of spee-ds are not considered in the. present^ paper . The 
purely supersonic case for small disturbances is also 
"tractable and flutter calculations for this case are being 
prepared. 


PROCEDURE 


In the idealized case of a wing with two degrees of 
freedom, wing bending and wing torsion, and based on two- 
dimensional air forces, the de terminantal equation yielding 
the flutter condition may be written' in the form (refer- 
ences 4_ an, l 5) 


R 

R 


aa 

-h 

^aa 

ca 

+■ 

ilea 


^ah^L ilah 
^ch + ilch. 


0 


The two real equations contained in the complex 
determinant may yield in any given problem the two 
unknowns , flutter speed and flutter frequency. 

Expressions far the R’s and I's in the incompressible 
case are listed following the definitions of the various 
symbols in the appendix, and the evaluation of the terms 
is facilitated by the use of table I . In the compressible 
case the. R's and I's are expressed in terms of - the nota- 
tion of Frazer and Skan (reference J ), and table 
tains values for M = 0.7. 
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Instead of a direct solution of the flutter speed 
and flutter frequency for a case in which the wing struc- 
tural parameters are given > it is more convenient to solve 
for the parameter ^h/^a {.the ratio, of frequency in 
bending to fr equency in torsion) whi cir be longs to the 
border-line case of flutter for a definite chosen value 
of the parameter l/k. The elimination of X from the two 
equations contained in the determinant yields a quadratic 

equation in (wh/^a)^ from which coh/co a - may be found 
and, subsequently, X may be evaluated; and X, together 
with the given value of l/k, determines the flutter speed 
and frequency. It is convenient first to perform the 
calculations for the compressible case with M -~0~.7 and 
finally, in order to furnish the desired numerical com- 
parisons, to perform the calculations for the incompressible 
case (M = 0) utilizing the given structural parameters 
and the derived values of the frequency ratio ^h/^a- 


RESULTS AND DISCUSSION 


The main numerical results are summarized in table III 
and are shown plotted in figures 1 to 1+.. (No tabular 
values are Included for fig. ii. ) The ordinate In the 
figures is the flutter-speed coefficient v/bco a where ba> a 
represents a convenient reference speed. The abscTssa'is 
the frequency ratio 

The parameter ic may be considered to determine the 
wing density at a given altitude; thus, sc = 0.025 rep- 
resents the highest wing density used and K = 0.2. the 
smallest. Alternatively, the change In sc may be Inter- 
preted to represent a change in altitude for a given wing, 
and a change from sc = 0.10 to K = 0.05 may be con- 
sidered to represent an. altitude change from sea level 
to an altitude at which p equals one half the density 
at sea level, or approximately 18,000 feet. The normal 
drop in sonic speed with altitude should be taken into 

account in interpreting M = X. (see reference , fig. 2). 

c 

The examples treated may be further classified by values 
of the parameter a representing^ positions of the tor- 
sional Elastic axis; thus, a = -6.I4., a = -0.2, and a = G 
represent, respectively, elastic axes at J>0 percent, 

I4.0 percent, and 50 percent chord from the leading edge. 
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Also xa = 0.2 represents a position of the center of 
gravity 10 percent of the chord hehind. the el&stxo axis. 

The figures show flutter curves calculated for the 
two vhlues (M- = 0 and M - 0.7) representing respec- 
tively a low-soesd or incompressible case and a high- 
speed or comnressible case. For the ushkl circumstance 
of low values of the frequency ratio £%/to a , the effect 
of compressibility on the' flutter speed is seen to be 
relatively small; for the lower wing densities, the effect 
is a small increase, whereas for the highest wing density- 
used the- effect is a small detrimental one. 

i^r the divergence speed (frequency a) — -£0) the 
formulas of the static case are applicable and the slope 
of ther lif t curve increases according to. the Glauert- 
Prandtl rule, which yields the approximate formula 



1 

(l - M 2 )^ 



for very heavy wings (k — > 0) the values of l/k 
for flutter appro-ach <»; that isr; - the low frequency or 

the static case is approached and the (l - M 2 ) ' rule 
given in reference if. appears applicable. The empirical., 
formula for the flutter speed of reference ii (p. 17 ), 
which -is valid for high wing density and low values 
of <0^7^, may be modified to read 


v f 

bw a 


(l - 



1/2 


+ £l + X 


a 


Since the. slope of the lift curve does not increase in 
accordance "with the Glauert-Prandtl formula beyond a 
certain value- of M < 1, the formula is clearly 
inapplicable beyond a certain value of M. This value 
may be taken roughly to be in' the order of & = 0.7 
to 0 .75. For M = 0 . 7 , the formula indicates a decrease 
in the flutter speed of approximately 17 percent. 
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The effects of internal damping and of the modes of 
vibration have been omitted in the calculations. Inclusion 
of these effects would tend to reduce further the differ- 
ences between the numerical results for the compressible 
and incompressible cases. This statement is borne out by 
the results of reference 5 for a tapered wing. 

Calculations for a wing with an aileron cannot be made 
by the method of reference 1 wit*hout very extensive and 
difficult computations. It may also be remarked .that the 
numerical tables of references 1 and 3 show the' need of 
additional extensions and recalculations. 


CONCLUDING REMARK 


The main conclusion to be derived from study of the 
numerical flutter calculations is that the effect of com- 
pressibility on the flutter 3peed (wing bending - wing 
torsion, no aileron) for subsonic speeds with no shocks, 
although complicated, is relatively small in the usual 
cases and, for a Mach number of 0.7 can be allowed for 
by corrections of small order to the incompressible-case 
results. . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 18, 1945 
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APPENDIX 

SYMBOLS AND NOTATION 

- ' ■ -«E 

half -chord used aa reference unit' length 

■ ‘ I 

coordinate of axis of rotation? (torsional axis) 
measured from midchord (reference 1|.) 

air density 

ratio of mass of cylinder of air of diameter equal 
to chord of wing to mass of wing, both taken 
for equal length along sp|n; this ratio may be 

expressed as ft = 0.2ij. wiiere ^ 

is weight - in pounds per foolnsp&n, b is in 
feety—and p/p 0 is ratio of actual air density 
to standard density at sea level 

location pf^ center of gravity of airfoil measured 
from a (reference l±) 

radius of gyration of airfoil referred to a 
(reference ij.) 

.... natural angular frequency of torsional vibrations 
around a in vacuum 

— natural angular frequency of wing in bending in 
vacuum 

speed of forward motion 

spaed of sound in undisturbed medium 

Mach number (v/c) 

flutter op critical speed 
.angular frequency of wing vibrations 
reduced frequency 

5 ‘functions of k (see references 1+ and 5 ) 



\ = 2k 
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In each of the following formulas the incompressible- 
case formula is listed first in the notation of references 4 
and 5 followed by the conversion of the formula to make 
use of the notation and numerical tables of reference 3 • 

(See tables I and II.) 
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where 
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The functions tabulated for convenience in tables I 
and II are the primed quantities given In the following 
definitions * 


Raa ~ Raa' + “ 


a 


^ah “ Rah’ " 


-a 


x 


Rea ~ Raa T “ 


a 


Fch “ Rch T + ^h^- “ k 


^aa “ ic I aa’ 


J ah _ k ^ah’ 


• ca 


i T , 

k x ca 


^ch 


k IctL ’ 


% 
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TABLE . I 


THE R ■ AND I ' EDUCTIONS FOR H = 0 


a 

11 

iA 

m 

H a a' 

R*h’ 

Ron’ 


I«a' 

I Ah’ 

loo' 

loh' 


10 


-0.37500 

-64.18650 

-0.50000 

-.74460 

-1.23300 

-1.72300 

169.34600 

I69.IOI40 

166.6120) 

160.12300 

1 

S 2.446 


H 

0 

-.16640 

-.49920 

-.83200 

-0.76200 

-.94840 

-1.28120 

-I.614/X) 

|> 1.66400 


5 

a 

-0. 37500 
-4.09274 
-11.47506 
-18.78650 

-0.50000 

-.58860 

-.76580 

-.94300 

37 .766OO 
37.67740 

37.50020 

37.32300 

a 

11 

-.72760 

0.56920 

-.15840 

1.45520 

H 

3.55353 


-O.37500 

-1.87035 

-4.84935 

-7.81275 

-0.50000 

-.51950 

-.55850 

-.59750 

15.47300 

15.45350 

15.41450 

15.57550 

m 

fl.00000 
J .89980 
<1 .77920 
[ .76500 

0 

-.13300 

■39900 

.66500 

1.13500 

1.00200 

.73600 

.47000 

1. 33000 


2.5 


-0.37500 

-1.14050 

-2.68200 

-4.23750 

-0.50000 

-.48250 

:fe 

8.13750 

8.15500 

8.19000 

8.22500 

} -°- m 

■ 


1.42500 

1.30000 

1.05000 

.80000 

■ 



a 

-O.57500 

-.B1757 

-1.72655 

-2.66730 

M 

4.88600 

4.92572 

5.00516 

5.08460 

mm 

11 

-■59790 

1.59500 

.99510 

1.19580 

B 


a 

-O.37500 

-•.41093 


°m 

1.01922 

1.17910 

”J> -0.79940 

H 

0 

-.10788 

-.32364 

-.53940 

1.87820 

1.7703? 

\m 

^J> 1.078B0- 


0.66667 

1 

-0.37500 ’ 
-.34014 
. .32454 
-.38108 

-0.50000 

-.40981 

-.22544 

.-.04907 

0.06123 

.15142 

-.33179 

.51217 

^J> -0.90187 

fl.00000 
J .81603 
> .51062 
L .28857 

0 

-.10420 

-.31260 

-.52100 

1 .94387 
1.63967 
1.63127 
1.42287 


4-0 

■ 

■ 

-0.37500 

-'llet 

-0.50000 

-.40577 

::K 

-0.18580 

-:B8 

•26535 

HI 

fl.00000 

J .81343 

L'Sil 

0 

-.10260 

-.30760 

-.51300 

1.96830 

1.66570 

1.66050 

1.45530 


m 



-0.37500 

-.30469 

-.22179 

-.21585 

-0.50000 

-.40370 

-.21135 

-.01892 

-0.29936 

-.20314 

-.01071 

.18170 

^j> -0.96216 

fl.00000 
J .81221 
S .49770 
[ .26457 

0 

-.10174 

-.30522 

-.50870 

1.97956 
1.67762 
i. (57434 

1.47066 

^j> 1.01740 
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TABLE II 

THE H» AID I' RJKCTI0K3 POH If = 0.7 


W 


* 


< 


A = 2k 

W 

a 

R s.a' 

Rkh’ 

H 6a' 

Boh' 

l»a' 

i«n' 

*ea' 

Ich' 

0.2 

, 10 

{? 

2.00100 

-17.97100 

-57.56700 

-96.68700 

-1.06200 

-1.65000 

-2.82600 

-4.ool)0o 

201.38200 
200 .79I1OO 
195.61600 
198 .44000 

1 

> 5.8B300 

r 1.79150 
) 2.20510 

1 3.11)870 
L4-24790 

0.02700 

-.16740 

-.55640 

-.94520 

-3 .96780 
-4.16220 
-4.55120 
-4.94000 

j> 1.94450 

o ,4 

5 ‘ 

{? 

0.89050 

-3.3055O 

-II.5575O 

-19.61850 

-0.95100 

-1.16750 

-l.64ooo 

-2.11300 

43.13200 

42.89550 

42.42300 

41.95000 

|> 2.36400 

f 1 -69465 
) 1.76805 
S 2.07025 
^2.48005 

0.05730 

-.10210 

-.42090 

-.73970 

-0.83100 

1-59400 

0.6 

5.33335 

l o- 2 

0.51)089 
-1.15911 
-4-49333 
-7 .73909 

-a. 86667 
-.97689 
-1.19733 
-1.41800 

17.97800 

17.B677B 

17,64711 

1742663 

1.10256 

fl.66314 
y 1.64420 
S 1.69300 
41-85713 

0.P8B93 

::5g» 

-.63207 

0.24407 

.09987 

-.16855 

-.47693 

\ 1.44200 

O.B 

2-5 

( ii 

l 

0.41231 

-49219 

-2.2699t 

4.00619 

-0.82537 

-.87725 

-.98IOO 

-1.08475 

9.92187 

9.87000 

9-76625 

9.66250 

|> 0.51887 

r 1.65403 
] 

0.12605 

-.01150 

-.28660 

-.56170 

0.72055 

.58300 

.30790 

.03280 

j> 1-37550 

1.0 

2.0 

l :■* 

0.40712 

-.15752 

-1.271)00 

-2.57320 

-0.80880 

-.83024 

-.87512 

-.91600 

6.47680 

1:® 

6.36960 

^j> 0.21436 

r 1 . 69^60 

O.I792O 

::gl 

0.96912 

.83276 

.56004 

.28752 

j> I.3636O. 

2.0 

1.0 

; ii 

l 

O.51805 

.38121 

.10213 

-.181)11 

-0.56^10 
- !5l836 

1.92266 

1.93160 

1.94950 

1.96740 

j> -0.08949 

Ji® 

S 1.23474 
L 1.15950 

0.5394a 

■39440 

.10446 

-.18560 

O.9968O 
.85180 
• 56180 
.27180 

\ 1.45000 

3.0 

0.66667 

r 11 . 

u 2 

-.08298 

-O.39IH . 
-.37364 
-.33872 
-30379 

0.85401 

.87148 

.96641 

■94135 

j> -0.17465 

S .87056 

L -780^7 

O.6528O 

•49073 

.19060 

-.11753 

1.03113 

.87706 

:lfe 

j> 1.54067 . 

k.Q 

— 

,0.5 

J‘ -°4 
l °" 

O.O808O 
.05917 
■ .00613 
-.05995 

-0.28700 

-.27670 

-.23869 

-.20549 

0,48699 

■50329 

.^6850 

|> -0.16302 

P1.0B573 
J .92521 
S .70333 
1 .61371 

‘m 

-I05EI5 

0.^870 

: 50280 
.17220 

j> 1.65300 

5.0 

0.1). 

J ii 

l 

0.12866 

.IO96O 

.06564 

.01389 

-0.14550 

-.13577 

-.11631 

-.09684 

0.32638 

■55557 

.37504 

j> -0.09733 

Pl'.oo88o 
J .84582 
S .62165 
l .53320 

“S s 

ll 2620 

0.82501 

.65557 

.51609 

-.02319 

> 1.6961)0 
J 
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TABLE III 
NUMERICAL RESUI/TS 
[r a 2 = 0.25j x a = 0.2] 


K 

H = 0.7 

M = 0 

iA 

V“a 

“/“a 


iA 

. u h/ u a 

“/“a 


• ' a . = - 0.4 

0.200 

< 

r 1.250 

2.000 
2.500 
3.533 
1 3.850 

1.000 

.590 

.436 

.234 

0 

1.080 
.736 
.623 
. -507 

.450 

1.350 

1.472 

1-559 

I.69I 

1.740 

0.85 

1.63 

2.05 

2.62 

2.99 

1.000 

:8? 

.234 

0 

1.16 
.82 
•72 
.62 
• 57 

SB 

0.100 

< 

- l.ltfJO 
2.000 
2.500 
3-333 
L 4.990 

1.000 

.653 

.519 

.3.61 

0 

1.070 

Ml 

1.500 
i: 1% 

1.847 

2.130 

1-79 

2.23 

2.85 

5.60 

0.655 

•519 

.561 

0 

0~86 

.56 

1.54 

1.71 

1.90 

2.16 

0.050 

< 

1.700 
2.000 
2.000 
2.500 
2.500 
3.333 
5.000 
_ 6-490 

1.1O0 

.818 

2.002 

.690 

6.106 

•494 

.292 

0 

1.150 
.901 
■ 1.970 
•757 
5-785 
.627 
•490 
• 4 io 

1-950 
1.803 
3.941 
1.893 
14.462 
2.091 
2.448 . 
2.630 

1.13 

1.77 

2.20 

B§ 

5.4.0 

1.100 

.818 

.690 

•494 

.292 

O' 

1.27 

1.00 

"89 

":Pi 

.54 

1.44 

1.77 

1.95 

2.28 

2.63 

2.90 

0.025 

\ 

2.500 
2.500 
3.333 
3.333 
5.000 
„ 9.860 

0.889 

2.270 

.435 

0 

0.943 

1.295 

.724 

2.120 

.560 

.350 

2.358 

3.257 

2 :. 4 l 2 

7.068 

2.799 

5.420 

2.30 

i.ei 

3.34 

4.80 

7.84 

0.&89 

•435 

0 

1.04 

■ 3 
~~66 
• 51 

1:11 

2.72 

3.19 

4.00 

a a -0.2 

0.200 

< 

1.030 
2.000 
2.000 
2.500 
2.500 
_ 2.980 

1.000 

;446 

f:ig 

4.791 

0 

i.050 
.622 
3.107 
, -515 

4 :& 5 

1.080 

1.243 

6.215 

1.288 

11.363 

1.310 

m 

1.000 

•446 

.269 

0 

1.17 
• 77 

”68 

.61 

0.62 

i-.i 4 

1.23 

1.28 

0.100 

< 

r 1.350 

2.000 
2.000 
. 2.500 
2.500 

lm 

1.000 
• 574 
2.365 
•435 
3.922 

.253 

0 

1.030 

.698 

2.251 

m 

.410 

1.390 

1.396 

9.122 

1.5.99 

1.650 

0.65 

1.62 

1.99 

2.49 

2.89 

1.000 

•574 

•435 

•253 

0 

H 

IBjl 

0.050 

< 

~ 1.700 
2.000 

2.000 
2.500 
2.500 
3.333 
3.333 

5.000 
^ 5.700 

1.000 

•755 

1.347 

I.’® 

An 

.204 

0 

1.000 

.828 

1.637 

2i$ 

.423 
• 570 

1.690 

1.657 

2.647 

1.737 

4.094 

1.890 

8.027 

2.115 

2.130 

1.22 

1.67 

1.25 

2.20 

2.86 

HI 

1.000 

-755 

1.347 

• 594 

• 451 
.204 

0 

1.09 
.9 6 
1.44 
.82 

‘T70 

”58 

•54 

1.33 

1.60 

1.79 

1.80 

2.00 

2,21 

2.29 

0.025 

< 

2.500 

2.500 

3.333 

3.353 

5.000 
5.000 
„ 8.430 

0.838 

1.050 

iS 

.368 

2.279 

0 

O.875 

1.047 

1I364 

.503 

2.010 

.540 

2.189 

2.619 

2.221 

4-547 

2-517 

10.048 

2.670 

B1M 

0.838 

1.050 

.588 

.368 

0 

O.97 

1.15 

•77 

.61 

•49 

2.20 

2.24 

2 - 49 

2.83 

3 - 17 
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TABLE m - Concluded! 
NUMERICAL RESULTS - Concluded! 
|tq 2 = 0.25; * a *0.2] 





M = 

0.1 



M = 

0 


‘ 

-iA 

V“a 

u>/cc a 

T/W 

1A 

idh/fea 

“/"a 

v/bdUa 






a = 0 



" 




' 1.150 

1.500 
2.000 
2.000 

2.500 
2.500 
3-333 
5.000 

„ 2.610. 

1.000 

.600 

lilt 

1.020 

.710 

liS 

1# 

1.500 

1.512 

.430 

1.170 

1.060 

1.086 

2.832 

1.110 

X.fS, 1 , 

0.12 

.85 

1.23 

1.000 

.600 

.343 

1.16 

.90 

•77 

0.14 

% 

0.200 

< 

.126 

1.518 

1.570 

1.590 

0 

1.45 

.126 

.70 

1.01 



5*000 

7.560 

1.110 













1.48 

0 

.69 

1.02 



' 1.350 

1.000 

.800 

i:SI- 

All 

l.Olff 

.850 

.637 

1-552 

i:81 

.L2L 

1.617 

1.764 

.380 

1.360 

— 


— 

— 

0.100 

< 

2.000 

2.000 

2.500 

2.500 

3.333 

3.333 

5.000 
„ 3.690 

1.275 

m 

f:gt 

3.111 

1.420 

1.43 

i-76 

0.515 

.373 

0.83 

"74 

1.18 

1.29 



■Ilk 

1.752 

1.943 

0 

2.17 

.174 

.64 

1.40 



2138 

0 

“Tio 

1.43 

0.050 

J 

' 2.000 
2.000 
2.500 
2.500 
3.353 

0.700 

1.092 

ill® 

.385 

1.475 

.147 

1.757 

0 

0.768 

1.070 

.642 

1.204 

■517 

1.361 

.373 

1.576 

.340 

1.556 

2.141 

1.605 

If 

1.68 

1.43 

2.15 

2.60 

2.74 

0.700 

1.092 

1 lie? 
.385 

0.90 

1.21 

.78 

1.51 

i:B 

f:? 9 x 



5.000 
5.000 
L 5.1*30 

1 M 1 

7.879 

1.870 

3.60 

3.91 

.147 

0 

■ 55 

•51 

1.96 

1.99 

0.025 

< 

■ 2.500 
2.500 
3-333 

3.535 

5.000 

5.000 

7.940 

0.797 

-38 

1.174 

xM 

0 

0.819 

.906 

.618 

1.101 

.458 

1.285 

.310 

m 

2.059 

i\m 

6.424 

2.450 

2.36 

2.31 

5.21 

4-44 

6757 

0.797 

:® 

.545 

0 

0.90 

.97 

.71 

.57 

’.42 

2.12 

i:i 

2.52 

2.76 
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Fig. 1 

























